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By expanding the wave function in plane waves, we use the pseudopotential method of density 
functional theory within the generalized gradient approximation to calculate the effective magnetic 
coupling energies of the S = 5/2 spins in the Fe2 dimer, approximated as Fe202H2. Setting the 
Fe-0 bond length at the value corresponding to the minimum total energy, we find the difference in 
antiferromagnetic and ferromagnetic exchange energies as a function of the Fe-O-Fe bond angle 6. 
The effective interaction is antiferromagnetic for 63° < 9 < 105°, and is ferromagnetic otherwise. 
Full potential augmented plane wave calculations were also performed at 6' = 100, 105°, confirming 
these results, and providing information relevant to the local anisotropy of the spin interactions. 

PACS numbers: 71.15.Mb, 61.46.+W, 75.75.+a 



Single molecule magnets (SMM's) have been under in- 
tense study recently, due to their potential uses in mag- 
netic storage and quantum computingJ^SiS, The mate- 
rials consist of insulating crystalline arrays of identical 
SMM's 1-3 nm in size, each containing two or more 
magnetic ions. Since the magnetic ions in each SMM 
are surrounded by non-magnetic ligands, the intermolec- 
ular magnetic interactions are usually negligible. Al- 
though the most commonly studied SMM's are the high- 
spin Mnl2 and FeSf^*^ such SMM's contain a variety 
of ferromagnetic (FM) and antiferromagnetic (AFM) in- 
tramolecular interactions, rendering unique fits to a va- 
riety of experiments difhcult^ 

In many single molecule magnets, and in Fe8 in par- 
ticular, the magnetic core contains [LM(OR)]^^ dimer 
ionS(^ where M denotes a magnetic ion (eg. Fe'^"'"), L is 
a ligand, and OR denotes an alkoxide ion, with R = H, 
CH3, CH2CH3, etc. In Fe8, the overall magnetic clus- 
ter is {[(tacn)6Fe802(OH)i2]Br7-H20}[Br-8H20], where 
tacn is 1,4,7-triazacyclononanei^ Near- neighbor Fe"^"*" 
ions are coupled via superexchange in four distinct pair 
bondings: four equivalent pairs of Fe'^"'" ions are coupled 
through two (OH)" ions with Fe-O-Fe angle 9 100.7°, 
one central pair is coupled via two O^^ ions with Fe- 
O-Fe angle 9 = 104.4°, and each of those Fe-^^ ions is 
coupled via the same 0^~ ions to two other Fe'^^ ions, 
and finally, four equivalent pairings via a single (OH)^ 
ion^ii Here we focus upon the magnetic superexchange 
interaction between the constituent magnetic ions that 
is mediated via two oxygen ions, which is different that 
the usual case of superexchange via a single oxygen ion. 
The attachment of alkyl R groups to the oxygen ions is 
probably of minor importance, because the oxygen or- 
bitals involved in the OR bond are orthogonal to those 
involved in the Fe-0 bonding. The local magnetic order 
depends strongly upon the Fe-O-Fe bond angle 9 and the 
bond length, which are determined by the local ligand 
environment. 

In addition, Le Gall et al. synthesized and measured 
the magnetization of four species of the isolated dimer 
Fe2,* and electron paramagnetic resonance (EPR) exper- 
iments were performed on one of them.— These dimers 



have the general formula [Fe(OR)(dK)2]2, where dK is 
a /3-diketonate ligand. In these dimers, the oxygen ion 
in each alkoxide group forms a bridge between the Fe"^"*" 
S=5/2 spins, allowing them to interact via superexchange 
through both 0^~ ions, as in the [LFe(0II)]2^ ions 
present in FeS. From magnetic susceptibility measure- 
ments, fits to the isotropic Heisenberg exchange model 
TL = — J Si ■ S2 were made.® Note that we use the sign 
convention in which a positive J corresponds to ferro- 
magnetic couplings. In comparing these four Fe2 dimers, 
they found no correlation between J and the average Fe- 
O bridging bond distance, but a linear correlation was 
found between J and the Fe-O-Fe bond angle 9, with 
J decreasing monotonically from -14.8±0.5 to -19.0±0.6 
cm~^ as 9 increased from 101.8° to 103.8°, respectively!^ 
One of these compounds, [Fe(OMe)(dpm)2]2, was studied 
with Mossbauer spectroscopy.- ° 

Le Gall et al. also compared these trends with predic- 
tions based upon extended Hiickel calculations performed 
on the simpler model system, [Fe(OH)H4]2r' using the 
approach of Hay, Thibeault, and Hoffman that for mag- 
netic dimers, J oc '^ii^Ei)'^ , where AEi is the energy 
separation between symmetric and asymmetric combina- 
tions of coupled magnetic orbitals^ Unfortunately, when 
they used the oxygen orbital parameters normally ex- 
pected, the J values calculated in this way increased with 
increasing 9, opposite to the experimental observations 1^ 
That approach was also not applicable for ferromagnetic 
exchange couplings. 

The static and dynamic properties of Fe2 dimers were 
studied theoretically by Efremov and Klemm, and inter- 
esting low temperature quantum steps in the magneti- 
zation were predicted. Recently, those authors studied 
local spin anisotropy effects, and found that the details of 
the low temperature quantum magnetization steps could 
be complicated by such local anisotropics.^^ 

Because of the very weak antiferromagnetic inter- 
actions (J « -0.9 ± O.lmeV) between the Fe(III) 
spins in the two related dimer compounds r^iiS4iSiiLi& ^- 
oxalatotetrakis(acetylacetonato)Fe2 and [Fe(salen)Gl]2, 
where salen is iV, iV'-ethylenebis(salicylideneiminato), 
magnetization steps were observed at low temperatures 



2 




FIG. 1: Sketch of the [HFeO]2 dimer. The Fe-O-Fe angle is 
denoted 6. 



in themii2i2£ In the former case, all five of the magnetiza- 
tion steps were observed using pulsed magnetic fields, 
and they were found to be evenly spaced, suggestive 
of an isotropic antiferromagnetic Heisenberg exchange 
interaction.^^ In the second case, the lower two and part 
of the third magnetization steps were observed by pulsed 
magnetic fields at low temperature. Although the sec- 
ond step was rather sharp, with a sharp dM/dH peak, 
the first step had a much broader lincwidth, suggestive 
of two primary constituentsi^ Although the sample mea- 
sured contained many unoriented crystallites, the broad 
first peak followed by the sharp second peak is consistent 
with a substantial amount of local spin anisotropy of the 
type -Ja ELi where J, « O.ljM. 

Here we focus on the constituent dimers present in 
Fe8, which are of the [RFeO]2 type, in which each Fe''+ 
ion shares electrons with one aliphatic R group and two 
bridging ions, which have no aliphatic substitutions 
on them, as pictured for R = H in Fig. 1. We model this 
system by replacing the R group with H. To study the 
geometry effect upon the magnetic interaction and the 
spin state of the dimer, we perform a first principles cal- 
culation of the electronic and magnetic structures and of 
the total energy of the [HFeO]2 dimers for ferromagnetic 
and antiferromagnetic states and for 57° < < 110°. 
From these results, we are able to obtain J{9). 

The majority of the calculations presented here were 
performed within the density functional theory, using 
the generalized gradient approximation (GGA) for the 
exchange-correlation functional, and the pseudopo- 
tential method combined with the expansion of the 
wave functions in plane waves. We assume a three- 
dimensional structure of such [IIFeO]2 dimers, and for 
simplicity, pack them into a periodic array of tetragonal 
unit supercells of dimension 12x12 x IsA. This unit su- 
percell is sufficiently large that the interaction between 
neighboring dimers is utterly negligible. For all atoms in 
the dimer, we use the ultrasoft pseudopotentials.^^ We 
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TABLE I: Magnetic characteristics of the [HFeO]2 dimer cal- 
culated using the pseudopotential method for different Fe-O- 
Fe angles 6. 

set the cutoff energies for the plane-wave expansion at 
400 eV, and perform the calculaton only for the F-point 
in the first Brillouin zone, which is sufficient for a single 
molecule calculation. To obtain the equilibrium struc- 
ture of the dimer, we apply the optimization procedure 
that relaxes the system until the forces acting on each 
atom converge to within 0.02 eV/A. We obtain the val- 
ues Sfm and Safm for the local spins on the Fe sites 
in both the ferromagnetic and antiferromagnetic config- 
urations by integrating the local spin densities over a 
sphere of radius 1.1 A centered about the Fe site. For 
each angle 9, the total energies EpM and Eafm foi' the 
antiferromagnetic and ferromagnetic configurations are 
calculated. Then, the exchange coupling J is determined 
from 

J _ Eafm — Efm ^-^-^ 
SafmSfm 

Ideally, if the absolute energy values EpM and Eafm 
were reliable, we would calculate J from Eafm I Safm ~ 
Efm / Spj^j. We make the above approximation because 
the absolute energy values are not nearly as reliable as 
are the much smaller energy differences. 

In order to find the equilibrium structure of the dimer, 
we first minimize the forces acting upon the atoms within 
the dimer. We find that the minimum energy is reached 
for antiferromagnetic states in which the Fe-0 bond 
length e = 1.83 A, and Fe-O-Fe bond angle 6 = 82.4°. 
We then keep the bond length fixed to that optimum 
value, and calculate Efm, Eafm, Sfm, and Safm for 
six 9 values in the range 57° < 9 < 110°. Our results are 
presented in Table 1. 

We find that for 9 > 105° and for 6* < 63°, the dimer 
prefers to be in a ferromagnetic ground state, whereas 
for 63° < 9 < 105°, an antiferromagnetic coupling is pre- 
ferred. The change in magnetic order is accompanied by 
strong changes in the spin density on the Fe ions, as well 
as in the spin density on the O ions that control the su- 
perexchange interactions, which data are not presented. 
For the Fe dimer components in Fe8, the [LFe(0H)]2^ 
effective dimers have 9 = 100.7°, which our calculations 
indicate is likely to be antiferromagnetic, contrary to the 
experiments. The central effective [LFeOJj^ dimer has 
9 = 104.4°, which is borderline between ferromagnetic 
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FIG. 2: Plot of the exchange energy J in meV versus the Fe- 
O-Fe angle in degrees. Shown are our results from Table 1 
(solid blue square), a guide to the eye (solid blue curve), and 
the data for Fe2 dimers (red crosses).— 



and antiferromagnetic in our calculations. Experimen- 
tally, it is ambiguous, because the ground state configu- 
ration of those spins are also affected by superexchange 
via a single oxygen in an 0H~ ion, which is likely to be 
strongly antiferromagnetic. 

In Fig. 2, we have plotted our J{9) data, along with 
the exerimental data of Le Gall et al. for the four Fe2 
dimers— We note that the experimental data points are 
close to the guide to the eye between our calculated angles 
= 100, 100°, so that our pseudopotential calculation is 
close to predicting the experimentally observed exchange 
energies of these four compounds. 

We remark that Pederson, Kortus, and Khanna have 
performed a much more accurate set of calculations for 
Fe8, and obtain the experimentally found S— 10 ground 
state, with the central [LFeOJj"'' dimer being antiferro- 
magnetically coupled to the other Fe spins^ They also 
found global anisotropy parameters in rather close agree- 
ment with experiment as deduced from magnetization 
experiments and from oscillations in the quantum tunnel 
sphttingiS*2& 

We note that although the dimer we have studied is 
rather different from the Fe2 dimers that exhibit superex- 
change via the oxygen ions in alkoxides, to the extent 
that the alkoxides can be approximated by the ions 
we studied, the 9 values observed in the four Fe2 dimers 
varied from 101° to 104°, which are in the antiferromag- 
netic regime.^ As in the extended Hiickel calculations, it 
appears that the general trend we calculate of increasing 
J with increasing 6 is opposite to the experiments on the 
four Fe2 dimers. However, we do not have enough data 
points between 100° and 105° to determine if the gen- 
eral trend is accurately followed. More important, there 
appear to be serious discrepancies with the interesting 
case of [Fe(salen)Cl]2, with 9 « 90° and a very weak, 
antiferromagnetic exchange interaction«iSiiS*iiiSi2S 



Fe d 


FM 


FM 


AFM 


AFM 


AFM 


AFM 


orbital 


Fel,2 T 


Fel,2 i 


Fel T 


Fel i 


Fe2 t 


Fe2 i 


d,2_,2 


0.8454 


0.2541 


0.8451 


0.2194 


0.2192 


0.8451 


dx2 — y2 


0.8822 


QAm 


0.8875 


0.0875 


0.0885 


0.8875 


dxy 


0.8953 


0.0549 


0.8930 


0.1321 


0.1312 


0.8931 


dxz 


0.7494 


0.2585 


0.8777 


0.3247 


0.3247 


0.8776 


dyz 


0.8910 


0.0562 


0.8923 


0.0921 


0.0918 


0.8923 



TABLE II: FLAPW results for the partial charges within the 
Fe d-orbitals for 9 = 100°. The FM and AFM cases for up 
and down spins within the Fel and Fe2 spheres are presented. 
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TABLE III: FLAPW results for the partial charges within the 
Fe d-orbitals for 6 = 105°. The FM and AFM cases for up 
and down spins within the Fel and Fe2 spheres are presented. 



In addition, our pseudopotential calculation does not 
contain any specific information regarding the relative 
spin densities in the various occupied Fe orbitals. Hence, 
we are unable to determine the local spin anisotropy 
parameters relevant for more detailed studies of Fe2 
dimerssi^ In order to obtain useful information regarding 
the local spin anisotropics, another calculational proce- 
dure, such as a full potential augmented plane wave cal- 
culation, would be necessary. In [Fe(salen)Cl]2, the un- 
usual behavior of the first magnetization step is strongly 
suggestive of a substantial local spin anisotropy interac- 
tion, as noted ahoveM^ 

As a preliminary check upon the validity of the pseu- 
dopotential method and also as a first-principles investi- 
gation of the possibility of local spin anisotropy, we have 
made full potential augmented plane wave (FLAPW) cal- 
culations of the spin densities and energies of the [HFeO]2 
dimer at = 100, 105°. In Tables II and III, we pre- 
sented our results obtained using the FLAPW method 
for the partial charges in the various Fe(III) d-orbitals 
for both up (f ) and down (|) spin configurations in the 
FM and AFM configurations, on Fel and Fe2 sites, for 
9 = 100, 105°, respectively. In Table IV, the energy dif- 
ferences between the FM and AFM cases for these angles 
are also presented. 

In comparing our FLAPW results from Table IV with 
those obtained from the pseudopotential method in Table 
I, we see that they both predict a crossover from AFM 
behavior for 9 < 100° to FM behavior at larger angles. 
Although we did not perform the pseudopotential cal- 
culation at 6* = 105°, our FLAPW calculation at 105° 
predicts FM behavior, so that the predicted crossover at 
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6 (deg) Eafm - Efm (eV) 
100 -0.69134 
105 1.0594 



TABLE IV: FLAPW calculations of the energy differences 
between AFM and FM states for 6 = 100, 105° 
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TABLE V: Total charge densities on the Fe 1 and 2 sites 
for up and down spins, respectively, for the FM and AFM 
configurations, at the two angles studied using the FLAPW 
technique. 



9c is for 100° <9c< 105°. In addition, it is evident that 
the largest differences in the spin densities at within the 
dxy and dyz orbitals for the FM case for both 9 values, 
whereas comparable (for = 100°) or slightly larger (for 
9 = 105°) differences in the spin densities occur within 
the (1^2 _y2 orbitals for the AFM configurations, respec- 
tively. We expect that the bonding orbitals will be the 
dxz and d^^-r^ orbitals, for which the magnitudes of the 
differences in the spin densities were found to be nearly 
equivalent for both the AFM and FM configurations, for 



both angles studied. Thus, the spin density differences 
within the bonding orbitals, that participate in the ex- 
change interactions, depend very weakly upon the Fe-0- 
Fe bond angle 9. The greatest 9 dependence is for spin 
density differences in the non-bonding orbitals, and the 
orbital that shows the greatest distinction between FM 
and AFM behavior is the d^2_j,2 orbital. Thus, it appears 
that the da;2_j,2 -orbital governs the local spin anisotropy 
in Fe2 dimers. 

In order to estimate the strength of the effective 
Hciscnberg coupling, if we take the differences in the to- 
tal spin densities from Table V, we obtain J ~ —57 mcV 
and J = 86 meV for 9 = 100, 105°, respectively. Taking 
the spins to have the value 4 in each case, we get J = — 43 
meV and J = 66 meV, respectively. These numbers are 
larger in magnitude that the values obtained using the 
pseudopotential technique. 

We remark that the most interesting Fe2 dimer to date 
is [Fe(salen)Cl]2, which has CI ions in place of the H 
ions pictured in Fig. 1. In addition, there are organic 
ligands attached to the O ions, as in all of the Fe2 dimers 
presently known&iSiSS We have not yet studied this very 
interesting case, but plan to do so soon. The Cl^ ions 
could polarize the local spins in the Fe c?-orbitals, and it 
would be interesting to see if the d^^-r^ orbitals would 
be the most affected. 
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